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Abstract 



We examine a simple heuristic test of integrability for quantum chains. This test is 
applied to a variety of systems, including a generic isotropic spin-1 model with nearest- 
neighbor interactions and a multiparameter family of spin-1/2 models generalizing the 
XYZ chain, with next-to-nearest neighbor interactions and bond alternation. Within the 
latter family we determine all the integrable models with an o(2) symmetry. 
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1. Introduction 

This work examines integr ability criteria for quantum chains. In particular, we pro- 
pose a simple heuristic test of integrability for quantum chains with short-range interac- 
tions, consisting essentially in proving or disproving the existence of a nontrivial conserved 
local charge involving three-point interactions. We indicate the conditions in which this 
might be a necessary condition of integrability, as well as circumstances in which it is 
expected to be a sufficient condition. The test is then applied to a number of models. Our 
approach is motivated by our recent extensive analysis of the structure of the conservation 
laws for integrable spin chains of the XYZ type and the Hubbard model [1]. 

For hamiltonian systems, the common definition of quantum integrability mimics the 
Liouville-Arnold definition of classical integrable hamiltonian systems: a quantum system 
with N degrees of freedom is called integrable, if it possesses nontrivial, functionally 
independent and mutually commuting conservation laws. 

For classical continuous systems, the proof of integrability usually amounts to display- 
ing a Lax or zero-curvature formulation. Although this is not always easy, there are various 
other manifestations of integrability that can be probed, such as a bi-hamiltonian formu- 
lation, nontrivial symmetries, prolongation structures or higher order conserved charges 
(see e.g. [2]). Furthermore, one can apply a systematical integrability test, based on the 
Painleve property [3]. For classical discrete systems, there are related integrability indica- 
tors [4]. But except for higher order conservation laws, these integrability signals are no 
longer available for quantum chains. 

On the other hand, the integrability of quantum chains is usually demonstrated rather 
indirectly, by showing that the model can be solved by the coordinate Bethe ansatz or that 
the hamiltonian can be derived from a commuting family of transfer matrices related to the 
Yang-Baxter equation. But these are only sufficient conditions for integrability. Moreover, 
testing these sufficient conditions is often not easyS 

It is thus clearly desirable to design a more general, simple and efficient integrability 
test for quantum chains. Here we propose a simple test based on the existence of a local 
nontrivial three-point charge Hs, the higher order conservation law just above H2, the 
defining hamiltonian of the model. Locality means that the interaction involving a certain 
set of sites disappears when the distances between them become sufficiently large. We will 

^ Note however, that the Yang-Baxter equation impHes a relatively simple equation, known as 
the Reshetikhin condition. See section 2 for a discussion of this point. 
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indicate below why for quantum chains, one can expect that the existence of should be a 
necessary and in some cases a sufficient condition for integrability. From the computational 
point of view, the advantage of this simple-minded approach is that the number of possible 
candidates for the three-point charge is usually not exceedingly large, and they are often 
restricted by the symmetries of the system. 

A related integrability test has been studied in [5]. However, these authors confined 
their test to the search of one nontrivial conservation law (not necessarily Hs ) in spin 1/2 
models, with the simultaneous existence of a ladder operator providing a recursive scheme 
for the calculation of the other conservation laws. This is certainly less general and less 
constructive than the test proposed here. 

Another integrability test has been formulated for a very special class of self-dual 
quantum chains. For such systems, there exists a simple sufficient condition for integrabil- 
ity, due to Dolan and Grady [6]. This condition actually applies to any type of self-dual 
systems, discrete or continuous, and defined in any number of space-time dimensions. 
These systems are described by an hamiltonian of the form: 

H2 = aH + PH; (1.1) 

where H is the dual of H, with duality being defined as any nontrivial linear operator with 
the property H = H. If such hamiltonians satisfy the relations: 

[H,[H,[H,H]]] = IQ[H,H], (1.2) 

there exists a (potentially infinite, for infinite systems) family of conservation laws and 
these can be constructed systematically. We stress that this sufficient condition for inte- 
grability is simply the requirement of the existence of a charge H^, of a particular form. 
We have thus a neat situation here in which the the existence of a third-order charge ifs 
guarantees the existence of an infinite family of conserved charges. Unfortunately, the 
applicability of this condition seems rather limited.@ For example, the Hxyz hamiltonian 
is self-dual, as it may be put in the form: 

HxYZ = aHYz + f3HYz; (1-3) 

^ The condition (1.2) appears to apply only to the Zn generalization of the Ising model defined 
in [7]. Systems satisfying (1.2) are sometimes called superintegrable and the underlying algebraic 
structure is the so-called Onsager algebra [8]. 
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where 

Hyz = Et^^^M+i + ^-^l^l+i]' (1-4) 

where cr" are PauH spin matrices acting nontriviaUy only on the site i of the lattice A. The 
duality is defined by: 

However, the condition (1.2) is not satisfied unless Xy = (in which case (1.1) reduces to 
the Ising model hamiltonian), or Az = (the XY model). 

Our proposal is described in detail in the next section. It is formulated for the class 
of translationally invariant models with nearest-neighbor interactions. These models are 
defined on a lattice A which is either finite with periodic boundary conditions (i.e. A = 
{1,...,A'"}, with A'" + 1 = 1) or infinite (i.e. A = Z). Notice that the restriction to 
nearest- neighbor interactions is not as severe as it might appear at first sight, since this 
class of models actually contains, when appropriately reformulated, any model involving 
binary interactions with finite range. Indeed, any model with a finite interaction range ko 
(where /cq = 2 corresponds to nearest-neighbor interactions) can be equivalently described 
in terms of nearest-neighbor interactions just by grouping together ko consecutive sites 
into a single one on which a vectorial spin-like variable would live. Similarly, any model 
with a more general invariance under a shift j — > j +jo can be made invariant with respect 
to a translation by a single unit, by grouping together jo consecutive sites. 

Let {S"'} denote a set of quantum operators (distinguished by the index a), acting 
nontrivialy only in a Hilbert space Vi. (We will assume in this work that all the Vj's are 
finite-dimensional; in general they don't have to be identical). The full space of states of 
such chain is then the tensor product (^j^y^ ^3 • '^^^ class of hamiltonians under consider- 
ation has then the following form: 

H2 = J2i93 + hj,j+i], (1.6) 

where gj = g{Sj) describes interactions at site i and hjj+i is some site-independent func- 
tion of Sj,Sj+i which describes nearest-neighbor interactions. 



2. Conjectured integrability tests for quantum chains 

2.1. A conjectured necessary condition for integrability of a quantum chain 

From the early days of sohton theory, it has been clear that the existence of a non- 
trivial conservation law (beyond those associated with the standard conservation of mass 
or charge, momentum and energy) is a very strong indication of the existence of an infi- 
nite number of additional nontrivial conservation laws. But for continuous systems, either 
classical or quantum, one cannot focus the attention on the existence or non-existence of 
the conservation law with a degree just above the one which usually plays the role of the 
hamiltonian. The reason is simply that symmetry considerations can prevent the existence 
of Hn for a set of values of n. Take for instance the KdV equation, 

Ut + Uj;j:x + 6uu^ = 0, (2.1) 

for which u has degree 2 in the normalization where dx has degree 1. Due to a hidden Z2 
invariance, there are no conservation laws with densities of odd degree. More generally, for 
Toda systems related to Lie algebras, the values of n for which is non-zero are related to 
the exponents of the corresponding Lie algebra [9] . But notice that when the conservation 
laws are calculated by a recursive method, the conserved densities associated to charges 
that vanish by symmetry are not found to be zero; rather they are total derivatives. 

As indicated above, these considerations apply for both classical and quantum systems. 
But in the quantum case, there is another possible source for the absence of a conserved 
charge of given order, which is that the density can be exactly proportional to a null vector. 
Again this is observed in the quantum KdV case, where for special values of the central 
charge, the free parameter in the defining commutation relation, some conserved densities 
become proportional to null vectors of the Virasoro algebra [10]. 

Now, our point here is that both these possibilities are absent in the case of quantum 
chains. In particular, it is clear that on the lattice there is no room for total derivatives. 
Furthermore, null vectors appear not to be relevant if the spaces Vi are finite-dimensional.l 

^ Constraints imposed on the factor spaces Vi (e.g. realized via a projection onto some 
subspace of Vi) may "project out" some charges of an integrable model - in other words some 
of the conserved charges evaluated in a restricted space of states may conceivably vanish. Such 
constraints would thus have an effect similar as null-vectors. However, situations like this can be 
avoided by considering the conservation laws in the full (unrestricted) Hilbert space. 
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On the other hand, for systems whose integrabihty can be traced back to commuting 
transfer matrices, the conserved charges are obtained from the expansion of the logarithm 
of the transfer matrix [11] and is never absent. This can be easily seen for fundamental 
systems [12], characterized by a Lax operator L„(A) proportional to the model's R matrix 
Rno{X), normalized so that i?no(0) = Pno- 

Lr,{X) = RnoiX). (2.2) 

Here we use the standard notation in which Rij acts in the tensor product space Vi ®Vj, 
where Vi is the vector space associated to the site i. The index zero refers to the quantum 
space Vo on which the matrix entries of Lj act. Pno denotes the permutation operator. 
(For simplicity all the vector spaces Vi are assumed to be equivalent here). The transfer 
matrix is defined as: 

T(A) = Try„L^(A)...Li(A). (2.3) 

For the fundamental models, a simple calculation shows that the second logarithmic deriva- 
tive of the transfer matrix at A = is proportional to: 

Qs ~ ^[hjj+i, hj+i,j+2\ (2.4) 

where is defined in (1.6). (Here we follow the notation used in [1]: Qn stands for the 

(n — l)-th derivative of the logarithm of the transfer matrix and is obtained from by 
stripping off contributions of the lower order charges.) If one excludes a pathological case 
in which all adjacent links commute,^ the three-point charge defined above is nontrivial. 
The existence of a nontrivial third order charge can also be proven for a more general class 
of integrable models, characterized by the condition 

L^{0) = PnO- (2.5) 

Note that such L matrix is not necessarily itself a solution of the Yang-Baxter equation. An 
example of a transfer matrix which satisfies (2.5), but which does not define a fundamental 

For example, such pathological situation arises for the "chopped XXZ" model introduced in 
section 3. 
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model, is provided by the transfer matrix of the Hubbard model found by Shastry [13]. For 
models satisfying (2.5), the second logarithmic derivative of the transfer matrix becomes: 

jeA jeA 

^ d^L ^2.6) 

+ J^T-' (O)Trvo (Pno . . . ^(0) . . . Pio)- 

The last two terms involve only nearest-neighbor interactions, while the first is a sum 
over triples of consecutive spins. This term has the same form as (2.4), and it is clearly 
nontrivial except again for pathological systems with all adjacent links commuting. 

The above considerations motivate the following conjecture for the class of systems 
described by (1.6): 

Conjecture 1: A translationally invariant periodic quantum spin chain, with a hamilto- 
nian H2 involving at most nearest-neighbor interactions, is integrable only if there exists 
a nonvanishing local independent charge H3, being a sum of terms coupling spins at most 
three sites, which commutes with H2 for all chain sizes N > 3. 

This conjecture implies then a simple test consisting in establishing the existence of 
such H^. In particular, one may conclude that a spin chain is nonintegrable, by demon- 
strating the non-existence of a nontrivial H^. On the other hand, if such a charge exists, 
the system is likely to be integrable; but of course its integrability has to be proven in- 
dependently. Actually, it appears also that for a large class of systems (and in particular 
for all the models considered in sections 3 and 4), the mere existence of H3 is enough to 
guarantee the existence of an infinite family of conserved charges in involution. 

2.2. Clarifying comments related to the formulation of conjecture 1 

Some elements entering in the formulation of the conjecture deserve clarification. 

i) Independence of the charges 

The charge in the above conjecture should be independent of the hamiltonian and of 
possible charges of lower order (such as, e.g. the components of the total spin). For general 
quantum integrable models, the issue of functional independence can be quite complicated; 
it is indeed a major difficulty in formulating a general definition of quantum integrability 
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in a completely rigorous way.El However, functional independence of the charges can be 
usually easily verified for spin chains with short range interactions: the leading term of 
contains n adjacent interacting spins; the leading terms of and Hm for m ^ n being 
clearly distinct, the corresponding charges are linearly independent. Furthermore, such a 
cluster of n adjacent spins cannot be obtained from a product of lower order charges (since 
that would also generate terms with n non-adjacent spins), 
ii) Stability of the charge under a variation of the chain lenght 

The last item in the above conjecture ensures that the existence of should not be 
affected by a change N N + k, where k is an arbitrary positive integer. This stability 
requirement can be simply illustrated for the XXX model defined by the hamiltonian 

N 

H = ^(T,-cT,+i (2.7) 

i=l 

(where as usual periodic boundary conditions are assumed). This hamiltonian commutes 
with any component of the total spin 

N 

5" = ^<. (2.8) 

1=1 

A simple calculation shows that for = 4 the quantity 

N 

H' = Y^a,-a,+2 (2.9) 

1=1 

is conserved. H' is also conserved for N = 5 but not for any higher value of N. The reason 
for this behavior is that for N = 4 and 5, H' can be regarded as a nonlocal charge. Such 
nonlocal charges typically can be obtained from powers of the local ones. More exactly, 
for A = 4 we have the following identity (modulo an additive constant): 

S^S" = 2H + H'. (2.10) 

^ Dropping the requirement of functional independence does not lead to a meaningful definition 
of integrability [14]. In the context of spin chains this can be easily seen. Consider an isotropic 
chain, for which the hamiltonian commutes with all the components of the total spin. Arbitrary 
powers of any of the spin components yield then a set of mutually commuting conserved charges. 
Removing the requirement of functional independence in the definition of integrability would 
therefore render any isotropic spin chain automatically integrable. 



7 



Similarly, for A?" = 5 

S'^S'' = 2H + 2H'. (2.11) 

However, for AT > 5 H' is not related to 5°'<S°' or other nonlocal charges and it is no longer 
conserved. For A?" = 6 for instance, we have instead 

S'^S'' = 2H + 2H", (2.12) 

where 

TV 

H" = ■ + ICT^ ■ CTi+3] (2.13) 

i=l 

contains contributions with one and two holes. Thus the conservation of H' is not preserved 
under the change N N + 2; the same is true for H" . This is actually typical for all 
such "accidental" charges, whose conservation for some particular values of is due to 
an accidental identity (relating them to nonlocal charges), which is true only for these 
particular values of A'", 
iii) Locality vs nonlocality 

The above example illustrate another issue: it is not always easy to distinguish the nonlocal 
charges (e.g. H') from the local ones (e.g. H) for finite chains. For an infinite chain, 
local and nonlocal expressions are quite distinct (the latter contain interactions between 
arbitrarily distant spins). On the other hand, for the XXX model with A?" = 4, 5, one hole 

in the expression of a 2-spin conserved law reflects "nonlocality" ! However, as exemplifled 
above, the form of those nonlocal charges that can be written as powers of local ones is 
strongly A^-dependent; this property makes them easily detected. Notice also that the 
local charges for a finite chain may be defined non- ambiguously from the densities of the 
first N charges of the infinite chain. 

2.3. A conjectured sufficient condition for integrability based on the boost operator 

The integrable points in a multi-parameter space of general spin chain hamiltonians 
can usually be simply characterized by the occurrence of a dynamical symmetry. Related 
to such symmetry is the existence of a ladder operator B, acting on the conservation laws 
as 

[B,Hn]=Hn+l, (2.14) 
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where denotes a charge with at most n adjacent interacting spins B This motivates 
the formulation of a simple conjectured sufficient condition for integrability, based on the 
existence of a ladder operator B, in conjunction with the presence of a nontrivial charge 
H3: 

Conjecture 2: A translationally invariant periodic quantum chain, with a hamiltonian 
H2 involving at most nearest-neighbor interactions, is integrable if there exists an operator 
B such that for all chains of length > 3, [B, H2] is nontrivial and 

[[B,H2],H2]=0. (2.15) 

In contrast with the first conjecture, the approach here is more constructive: it in- 
dicates how can be built, i.e. as [B,H2]. This constructive aspect presupposes that 
we can easily guess the form of S. In all the cases we have considered, such a B turns 
out to be proportional to the first moment of an appropriately symmetrized form of the 
hamiltonian H2, i.e. with 

H2 = ^hjj+i, for hjj+i symmetric in j, j + 1, (2.16) 

jeA 

B is found to be 

5 = (2.17) 
The condition of commutativity of H2 and assumes then a particularly simple form: 

+ hj+i,j+2, [hj,j+i, hj+ij+2]] = 0. (2.18) 

jeA 

The derivation of this result is very simple. Starting from 

H3 = [B, H2] =^U - l)[/ij-i,j, hjj+i] + (j + l)[hj+ij+2. hj^j+i] 



jeA 

= - X] [hj,j+i^hj+i,j+2] 



(2.19) 



jeA 

one then enforces 



J (2.20) 

+ [[hj,j+i^ hj+i,j+2], hjj+i] + [[hj+ij+2, hj+2,j+s]:hjj+i] = 0. 



^ Note that we allow for the possibility of a linear combination of lower order charges Hm<j 
on the rhs of (2.14). 
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Using the Jacobi identity, the first term on the rhs can be written in the form 
— [hj-ij,hjj+i\,hjj+i\; with the shift j ^ j + 2, it cancels exactly the fourth term. 
Then, by shifting j by one unit in the second term, we recover (2.18). 

From the Jacobi identity, this condition automatically ensures the existence of a second 
nontrivial conservation law which commutes with H2: 

[m,H2] = [[B,H3],H2] = -[[H3,H2],B] - [[i^s, S], i^g] = 0. (2.21) 

However, showing that [H4^,Hs] = and that higher charges commute with H2 requires 
additional information. For example, if it is known beforehand that the commutant of H2 
is abelian, (2.18) actually implies the existence of an infinite tower of charges in involution. 

2.4- Relation between conjecture 2 and the Reshetikhin condition 

It is known (see [15]) that the existence of such a ladder operator is a direct conse- 
quence of the Yang-Baxter equation for nearest-neighbor interacting chains for which the 
transfer matrix is a product of R matrices (the so-called fundamental spin chains) . But the 
conjecture is a priori independent of the Yang-Baxter equation and in principle there could 
exists models satisfying (2.18) and not the Yang-Baxter equation. Furthermore, (2.18) is 
easier to test that the Yang-Baxter equation for the related R matrix. 

Actually, for fundamental spin systems, (2.18) can be viewed as a condition for the 
matrix 

i?(A) = P[I + XH2 + 0{X^)] (2.22) 
to be a solution of the Yang-Baxter equation. This approach leads to the condition 

[hjj+i + hj+ij+2, [hjj+i,hj+ij+2]] = Xjj+i - Xj+i j+2 (2.23) 

for some quantity X . This relation first appeared in [16] (see eq. (3.20)) and it is attributed 
to Reshetikhin. An explicit derivation can be found in [17]. The Reshetikhin condition 
(2.23) is nothing but the local version of (2.18). It appears to be an anticipation of the 
boost construction of conservation laws. It is also pointed out in [16] that this equation is 
not satisfied by all integrable systems (which is by now understood from the fact that not 
all such systems are fundamental) ; in particular this is the case for the Hubbard model (in 
agreement with the conclusion in [1] concerning the non-existence of a boost operator). 

In the rest of this work we examine the existence of a third-order charge for a number 
of models. Some of the calculations were performed using "Mathematica" . 
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3. Example 1: spin-1/2 chains with next-to-nearest neighbors interactions and 
bond alternation 

3.1. Definition of the model 

We consider a 10-parameter family of spin-1/2 models, which contain, in addition to 
XYZ-type interactions, next-to-nearest neighbor interactions, bond alternation terms and 
a magnetic field coupling term: 

H2 = J2{ E [Aa^K+i + A'a(-l)V«a^i + A'Xa^2]+K}- (3-1) 

jeA a=x,y,z 

This is the most general spin-1/2 model with interaction range shorter then two lattice 
units admitting bond alternations. To ensure translational invariance, we assume that 
the number of sites is even. These models can be equivalently represented as nearest- 
neighbor interactions on a lattice A', whose bonds correspond to nonvanishing interactions 
(see Fig. 1). 



2 k .„ 2 k +2 2 k +4 




Fig. 1. The lattice A' corresponding to the hamiltonian (3.1), whose bonds correspond 
to nonvanishing interactions. 

A lattice with such a "railroad trestle" topology has been considered in [18] (in the 
case where all the couplings are equal). Notice that the lattice in Fig. 1 corresponds also 
to a generalization of (3.1) admitting bond alternation for next-to-nearest neighbors. A 
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particular case of such bond alternation yields then the "sawtooth" topology, which has 
been shown to posses an exact valence-bond ground state [19]. 

The models in (3.1) can be equivalently described by a hamiltonian of the form (1.6). 
The structure of (1.6) is recovered if we express (3.1) in terms of the variables: 

S^^ or S,^( ]. (3.2) 

The family (3.1) contains many interesting systems, including some that are well 
known integrable models and some which are "exactly solvable" in some sense. In par- 
ticular, among the class of isotropic (globally su(2)-invariant) models satisfying h — ^ 
and 

for all o, the following special cases are covered by (3.1): 

(i) the Heisenberg (XXX) model (A' = A" = 0): 

H2 = E^^>^+i' (3.4) 

(ii) the staggered XXX model (A" = 0): 

H,^J2^X + X'{-iy]a^a-+„ (3.5) 
jeA 

(iii) the alternating XXX model (A" = A = 0): 

H, = 5^A(-l)^-^>«+i, (3.6) 

jeA 

(iv) the Majumdar-Ghosh model [20]: (A' = 0, A" = | A): 

j&A 

As is well known this model possesses an exact valence-bond ground state. 
For models invariant under global spin rotation around the ^-axis: 



coscK smo! \ I (Ji 



ay J \^ — sin CK cos a J ya^ ' ' 



(3.8) 
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(in which case the global su{2) invariance is broken down to o(2)), the coupling constants 
satisfy: 

Some interesting o(2)-symmetric models that can be obtained from a specialization of (3.1) 
are: 

(v) the XXZ model (X'a = \" a = 0): 

H2 = Et^-K^J+i + ^'^J+i) + ^-^l^l+i]' (3-10) 

jeA 

(vi) the Lieb-Schultz-Matis model with alternating Heisenberg and Ising bonds [21]: 

X"a = 0, = Xy = X'x = X'y = -, 

A, = (l + t/)/2, = (l-t/)/2, 

with the hamiltonian 

H2 = ^W2j(^2j+l + U(T^j+i(T2j+2\^ (3-12) 

(vii) the Hubbard model: 

Xx — Xy = X' X = X'y = X" z = 0, X"x = X"y = 1, 

X', = ±x, = u, 

with the hamiltonian 



(3.13) 



H2 = E[4^7+2 + ^74+2 + U[l + (-l)1cT|a|+i]. (3.14) 

The equivalence of (3.14) (for a lattice with an even number of sites) with the usual 
formulation of the Hubbard model 

jeA 

(where Sj and U are two independent sets of Pauli matrices at site i), can be seen by 
redefining the spin variables in (3.14) as 

4 ^^2.' (3.16) 
13 



3.2. Strategy of the test 

For the completely general anisotropic natural candidate for the third order 

charge involves only nearest neighbor interactions on the the lattice in Fig. 1 and has the 
form: 

jeA ai=x,y,z 

+ (41^3 + (-l)^«lla3)^?^;^2<^3 (3.17) 
+ {a^'\ia,a, + (-l)^«Ha3)<^^;^l^;^3 

where aala2as and aa%.^a3 are arbitrary coefficients. We search for integrable systems by 
imposing the condition of commutativity of with the hamiltonian. This leads to an 
over-determined system of equations for the set of parameters in H^. This system has 
nontrivial solutions only for special values of the parameters of the hamiltonian. However, 
the analysis of this system is very cumbersome, since in absence of additional symmetries, 
( p.lTp contains 216 free parameters! Henceforth, we consider only three special situations, 
in which the analysis simplifies significantly: the isotropic case, the o(2)-symmetric case, 
and the anisotropic case without bond alternation nor next-to-nearest neighbor interaction 
- that is, the XYZh model. 

One might consider other candidate charges, involving other triples of sites than those 
in (13.171). However, suppose that one cluster different from those appearing in (|3.17D is 



introduced. Then (considering an infinite chain), to cancel the new terms arising in the 
commutator with the hamiltonian, an infinite sequence of other clusters would have to be 
added, with the distance between spins in such clusters growing arbitrarily. This would 
violate the requirement of locality. (The above reasoning is not so obvious however when 
some of the couplings vanish.) 

3.3. The XYZh model 

We first consider (3.1) in absence of bond alternation and next-to-nearest neighbor 
terms. The most general candidate charge ifs, coupling three nearest- neighbor spins, is: 

^3 = E E «aia2a3^?<^ia;^2, (3.18) 
jeA ai=x,y,z 

which contains 27 free parameters. By enforcing the commutativity of this candidate Hs 
with the hamiltonian, we found that if /i 7^ and no two couplings are equal, the only 
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solution is: aa^ajag = for all triples 010203. Thus there is no nontrivial charge H^, for 
the anisotropic XYZ chain in a nonzero magnetic field. This suggests thus that the XYZh 
model is nonintegrable, in agreement with the fact that the Bethe Ansatz solution for the 
XYZ model cannot be generalized to the case with a nontrivial magnetic field. 

3.4- The isotropic case 

In the isotropic case, the most general third order charge involving only the nearest 
neighbors on the lattice in Fig. 1 has the following form: 

J 6 A 

(7 + (-l)^7')<^-+2^-+4], 

where a, /3, /3, 7 and their primed variants are arbitrary parameters. Solving the commu- 
tativity condition [H2, H^] = 0, we found that (for N > 8) nontrivial solutions (for which 
not all parameters of are equal to zero) exist only in two cases: for the XXX model 
(A' = A"), with the solution a 7^ 0, /3 = = /3 = /3' = 7 = 7' = 0, i.e. 

Hs = J2^eabc(r^a''^+ia;i+^, (3.20) 
jeA 

and for A = A', which corresponds to two decoupled XXX models on the even and odd 
sublattices (the solution is then 7,7' 7^ and all other parameters of equal to zero). 

These results suggest that the only integrable isotropic models within the family (3.1) 
are of the XXX-type. In particular, the Maiumdar-Ghosh model, the alternating XXX 
model, and the staggered XXX model@ all fail the test of the existence of and thus 
seem to be nonintegrable. 



^ Note that in the continuous limit the alternating term corresponds to tr g in the WZNW 
model [22] with a level one afRne su(2) spectrum generating algebra, g stands for the basic field in 
the WZNW model, a 2 x 2 matrix in the su{2) case, tr g turns out to be an integrable perturbation 
of this WZNW model [23] . However, the staggered XXX model in the continuous limit gives the 
WZNW model perturbed by both tr g and the marginal current-current term. Together, these 
two perturbations are incompatible with integrability. Note also that for A = ±A' the staggered 
XXX model degenerates into a pathological "chopped XXX" chain consisting of disjoint bonds 
(see section 3.5 for a discussion of a similar case). 
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For N < 8 there exist nontrivial solutions for of the form (|3.19 ). Let us iUustrate 



this in the case A^;' = 0, when the hamiltonian (3.1) contains only nearest and next-to- 
nearest neighbors, 

H2 = H + ^H', (3.21) 

where n is an arbitrary parameter, and H and H' are given by (2.7) and (2.9) respectively. 
For = 5 we find that the XXX charge (3.20) commutes with (3.21). This is a consequence 
of the "accidental" identity (2.11) holding for N = 5. In other words, for N = 5 the next- 
to-nearest interaction is a "nonlocal" charge, (related to the square of the total spin) , which 
commutes with all the XXX charges. For N > 5 this solution disappears (the three-spin 
XXX charge no longer commutes with the hamiltonian (3.21)). Similarly, for = 6 there 
is a one parameter family of solutions 

Hs = ifs^^^ + i^Fs,! - 1/3 ifi + iiu- z/)F3%, (3.22) 

where u is an arbitrary parameter, and! 

^3,1 = E^"''[^>'+2^ -+3 + ^>'+l^,V3], (3.23) 
^3,2 = E^"''<^?+2^-+4- (3.24) 

jeA 

For N = 7 there is another solution: 

Hs = ifs™ + fiFs,! - ^2F|,2. (3.25) 
Similar accidental nonlocal three-spin charges, whose form changes with exist in fact for 



all A^, but it is only for N < 8 that they "look local", i.e. can be put in the form (|3.19|) . 



3.5. The o{2) -invariant case 

We are searching again for a nonvanishing charge of the form ( |3.17| ) . The require- 
ment of the o(2) invariance imposes a number of restrictions on the parameters aa}a2a3 

ana. Olaia2a3 ■ 



'-yxy "^zxz '-^xyx ^zyz (^yyx (^xyy C^zzx 
Olxzz Olzzy OLyzz Otxxy OLyxx OLxxx ^yyy 



(3.26) 



F^ 2 is the symmetric part of the quantity ^3,2 introduced in [1]. 
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and 

Oiyxz Olxyzi Oizyx C^zxyi Oixzy Oiyzxj 

(3.27) 



C^yzy CtLxzxi Olzxx Olzyyi ^xxz Olyyz- 



As a result, the number of free parameters in is now decreased to 56. We then check 
whether there are nontrivial of this form commuting with the hamiltonian. For N = 8 
for example, this leads to a linear system of roughly 240 equations with 56 unknowns, which 
depends on the 6 free parameters of the hamiltonian (A^, A;^, A^', A2', A^", A2"). Again one 
is looking for the values of these parameters allowing for the existence of a nontrivial 
solution (persisting when is increased). This system is best analyzed with a computer. 
The results are given below.i 

If Xx" is not zero, a nontrivial solution for ifa exists only in two cases: 

(a) for two decoupled XXZ models on two disjoint (even and odd) sublattices: (A^^ = A^;' = 
Xz = A2'), 

(b) the Hubbard model: (A^ = A^' = Xz" = 0, A2 = ±Xz')- The solution for the charge H3 
is then 

Hs =l][2A2'(aJa|+iaJ+2 " ^M+i<+2) + A."(aJa|+2^J+4 - ^j^j+2^hd 
jeA 

+ Xz'il - (-)^KaJaJ+2^|+3 - aJaJ+2^|+3 (3.28) 
-a|aJ+iaJ+3 + a|aJ+iaJ+3)], 

which can be translated, using (3.16), into the usual expression for the third-order charge 
in the Hubbard model [13, 1]. 

If Xj' is zero there are more possibilities. Nontrivial solutions for exist in four 
cases. These are: 

(c) the XXZ model: A^' = Xz' = Xz" = 0, 

(d) the staggered XX model: A^ = Xz' = Xz" = 0, 

(e) the staggered XXZ model: A^; = Xz' = Xz" = 0, 

(f) the model with alternating XXZ and Ising bonds: A^; = iA^'. 

The three-spin charge obtained for the XXZ model is: 

ZT \ ^ r \ f „x „y „z „x „z I „z X „y ^z^y „x \ 

jeA (3.29) 



^ Note that trivial o(2)-symmetric models with hamiltonians involving only the z-components 
of the spin variables at each site are explicitly excluded from consideration. 
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in agreement with [11,1]. 

For the staggered XX model 

H, = 5^(A + A'(-l)^)(aJaJ+i + ajaj+j, (3.30) 

the three-spin charge obtained via our test is identical to the XX charge, that is 

^3 = 5^a(aJa|+iaJ+2 + ajaj+iaj^^) + /?(aJa|+iaJ+2 " (3.31) 

where a and 13 are arbitrary coefficients. Another indication of integrability of this model 
is provided by its continuum limit, which corresponds to a theory of free massive fermions 
[24]. Since we are not aware of an explicit proof of the integrability of the lattice model 
( ^■30| ) in the literature, we present a simple direct proof by exhibiting a family of mutually 
commuting conservation laws. These can all be expressed in terms of the densities: 

^n,j — ^3 + 1 • • • ^j+n-2^j+n-l' \6.6Z) 

defined for n > 2. In terms of these quantities, the scalar and pseudoscalar conserved 
charges [1] of the XX model are: 



hi+) — \^p^^ -I- pi' 



^yy 
jeA 



Zl^nj - e^^j n odd 
jeA 



(3.33) 



and 



We also define 



7 ( — ) \ ^ xy yx 

jeA 
jeA 



(3.34) 



(3.35) 



jeA 
jeA 

The conserved charges of the staggered XX model contain two families For n odd 

these charges coincide with the XX charges: 



-"2m+l ~"'2m+l' 
-"2m+l —"'2m+l 
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(3.36) 



For n even, the conserved charges for ( p.30| ) are: 



(— ) ( — ) / (— ) ( — ) ' (— ) (o.oYj 



Mutual commutativity of the charges -ffn^'* as well as their commutation with the staggered 
XX hamiltonian (for |A| even) can be verified directly as in [1]. Note also that the boost 
operator, 

B = ^J2^[Hcr]a^+i + «+i) + A'(-l)^(aJaJ+i + ajaj^,)], (3.38) 

has the ladder property: acting on ( |3.30| ) it produces the scalar part of (3.31). 
Under a spin rotation by 7r/2 around the z-axis restricted to odd sites, i.e.: 



2j + l ^ ^2j + l' "2j + l ^ "2j + l^ 



(3.39) 



the alternating part of the staggered XX hamiltonian transforms into the two-spin pseu- 
doscalar charge of the XX model, 

+ ^ (3.40) 

jeA 

where 

4"^ = E[^K+i-^M+i]- (3-41) 

jeA 

This is a special case of the Dzyaloshinski-Moriya interaction [25]. Notice that the 
transformation (3.39) can be interpreted as a duality in the sense of (1.1) if we define 
o^2j+i — '^^2j+ii ^2j+i ~ ~^^2j+i (where the factor i has been introduced in order to 
have a = a). Then = and ^ = ik2~^\ One may then consider a general 

hamiltonian: 

H = Ai4+^ + X2h'f^ + X3ki~^ + \Aki^\ (3.42) 

where Ai, A2, A3 and A4 are arbitrary constants. This hamiltonian is integrable, as can be 
seen from the existence of an infinite family of conservation laws, given again (for n odd) 
by (3.36). Notice that (3.42) is self-dual for A1A4 = A2A3; however, it does not satisfy the 
Dolan-Grady sufficient integrability condition (1.2). 
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For the staggered XXZ model, defined by the hamiltonian: 



H2 = $^[A'.(-l)^(aJcTj+i + a}aj+,) + X.a^a^^^, (3.43) 
jeA 

the three-spin charge obtained from the test has the following form: 

jeA (3.44) 
+ A.'(aJa|^,a^,-aJa|+,a;V2)]- 

This charge can be also obtained using the boost operator as [B,H2]. Note that the 
transformation (3.39) establishes the equivalence of the staggered XXZ chain with the 
model 

H2 = Et^M+i - + ^^^I^I+i]' (3-45) 

jeA 

which is a particular case of the XXZ model with Dzyaloshinski-Moriya interaction: 

jeA 

where Jx, Jz and D are arbitrary parameters. Integrability of ( |3.43| ) follows then from the 
integrability of (3.46), which has been proven in [26]. Note also that the model (3.46) is 
not an o(2)-invariant one, and by a spin rotation (3.8) with a suitably chosen angle a it 
may be transformed into the anisotropic Dzyaloshinski-Moriya system: 

H2 = Y,[D.cr]ay^, + Dyaya^^, + D^a^a^^^. (3.47) 
jeA 

The model with alternating XXZ and Ising bonds, 

H2 =5^[(A. + A/) ((72^^.(7^^.+! + a|^af^+i) + (A, + A,>|^.(7|^.+i 

jeA (3.48) 
+ (A, - A/) 

^2j + lC^2j+2 + -^^ ^1^1+2]' 



which is a slight generalization of the Lieb-Schultz-Mattis model, presents certain pecu- 
liarities. The model has been diagonalized in [21] (for X" ^ = 0, A^ = A'^ = ^,Xz = 
(1 -|- f/)/2, X' z = (1 — f/)/2). As observed in [21], a convenient basis is provided by the 
eigenstates of Lj = (T2j + cr2j+i. Consider then subspaces of the space of states corre- 
sponding to a particular sequence {Mj} of the eigenvalues of the third component of the 
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Lj's. Since the Lieb-Schultz-Mattis type hamiltonian (3.48) commutes with each of the 
Lj's it does not mix different subspaces; in other words it is block-diagonal in this basis, 
and it can be diagonalized separately in each subspace. The projection operators onto the 
subspaces corresponding to different sequences {Mj} provide then a mutually commuting 
set of operators, all commuting with the hamiltonian. Therefore, (3.48) satisfies the defi- 
nition of integrability given in the introduction; admittedly the nature of these conserved 
charges appears a bit unusual and somewhat trivial. 

The explicit form charge found via the the test is: 

Hs = J2[aL]alj+^alj+s + /^^l.-^li+i^I+i], (3-49) 

where a and P are arbitrary coefficients. Interestingly, not only does the above sum 
commute with (3.48), but each term in it is separately conserved. Clearly, all these terms 
can be expressed as linear combinations of the projection operators discussed above. Let 
pii^j) jgj^Q^g projection onto the states with Mj = rrij. Then c"|jCr|^_|_i = + 
Pj~^^ — Pj^^) aiid Lj — Pj^^ — Pj~^^ ■ In particular, in the 2-'^/^-dimensional subspace 
where all the Mj's are zero (which is the sector containing the vacuum), (3.49) vanishes. 

Clearly, the block-diagonal nature of the hamiltonian (and hence the existence of 
the set of commuting projections) is preserved by the addition to (3.48) of an arbitrary 
interaction involving only the ^-components of the spin variables. Finally, we note that 
a particular case of (3.48) with Xz — \J = \J' = provides a pathological "chopped 
XXZ" system, consisting of N/2 disjoint XXZ bonds, with all neighboring links trivially 
commuting. 

Summing up, it appears that the apart from the models (a)-(f) described above, 
all other o(2)-symmetric models within the family (3.1) are nonintegrable. Among the 
integrable models, there are three situations in which the first moment of the hamiltonian 
(the boost) acts as a ladder operator for conservation laws: the XXZ chain (cases (a) and 
(c)), the staggered XX model (d), and the staggered XXZ chain (e). 

4. Example 2: isotropic spin-1 chains 

Consider now a class of isotropic spin-1 chains with nearest neighbor interactions. The 
most general hamiltonian contains a bilinear and a biquadratic term: 

= Et-^"-^^! + /^('^^'^^i)']' (4-1) 
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where 5'"'s are the su{2) spin-1 matrices, acting nontriviaUy only on the j-th factor of 
the Hilbert space 

<S)jC^- For convenience we choose the representation in which is 

diagonal, i.e. 

/O 1 0\ . /O -1 

^^ = -=1 1 ,^^ = ^1 -1 1,^^=1 |. (4.2) 
V^VOIO/ v2\oi 



Using the identity: 



where 



(S^S-.f = - ^S^Sl (4.3) 



Df = S'^S] + (4.4) 



the hamiltonian can expressed as: 



= - f + ^DfDf^,]. (4.5) 



jeA 

The boost operator yields the following candidate for the charge: 

HsiP) = E e'^''[{-2 + 2(3- ^)S^S'^+,S]+, - ^DfS]^,Df^^ 



jeA 

•\ad T~\db nc i Qa T\bd T~\dc 



+ {-(3 + ^){DfDf_,,S'^^^^ + s';d]%,d%^)]. 



(4.6) 



The commutator [ii/'2(/3), -f^3(/3)] vanishes only for /3 = ±1 or /3 = cxd; these cases have 
been already identified as integrable in the literature. For [3 = —1, (|4.1| ) reduces to 
the isotropic version of the Fateev-Zamolodchikov chain, associated with the 19-vertex 
model [27], whose integrability follows from directly from the Yang-Baxter equation. For 
(3=1, (|4.1|) describes the Sutherland sw(3) symmetric chain [28], whose hamiltonian can 
be rewritten in terms of the Gell-Mann matrices t": 

^2(l)~E^.%%i- (4.7) 

jeA 

It can be solved by the nested Bethe Ansatz [28] and it is also related directly to the 
Yang-Baxter equation. Note that in this case, -^3(1) can also be written in the form [1] 

= E/"'%"^?+i^.V2- (4.8) 

jeA 
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In the limit /3 — oo, (O) reduces to a system with purely biquadratic interactions, whose 
integrability has previously been established in [28, 29]. 

For the spin-1 models ( [4.1| ) with finite /3 7^ ±1, the boost operator does not produce 
a conserved quantity. Furthermore, as we show below, there exists no nontrivial local 
conserved charge involving up to three nearest-neighbors. The most general expression for 
such a charge would have the form: 

TT \ ^\ „ ^abc na oh qc , „ abc r~iad r~idb qc i ^ ^abc oa r^bd T\dc 

+ ci7Dj^Sj^iSj^2 + (^8^jDj+iSj+2 + ci9SjSj_^_i + aio{SjSj^i)^]. 

where the a^'s are undetermined coefficients. Enforcing the commutativity of H3 with the 
hamiltonian, we obtain as usual a number of constraints on these coefficients. In particular, 
the vanishing of the four-spin terms in this commutator requires: 

a2 - a3 = = ae = a7 = as = 0, (4.10) 

and 

a2 = a4 = if (3 = 0, 

(4.11) 

ai = (I+4//32 -4//3)a4, a2 = {-1 + 2/ (3)a^ if (3=^0. 

The above conditions means that in order for H3 to commute with H2, its three-spin part 
must be proportional to the commutator of the boost and H2 (both for /3 = and /3 7^ 0). 
The vanishing of the terms with two and three spins in the commutator imposes two further 
restrictions. First, the two-spin part of Hs must be proportional to the hamiltonian H2. 
Second, unless = 1 or /3 = 00 the three-spin part of must be trivial. 

The nonexistence of for ( f4.1| ) with finite /3 7^ ±1, suggests that all these models are 
nonintegrable. In particular, the bilinear system (/3 = 0), as well as the one with /3 = 1/3, 
for which there exists an exact valence-bond ground state [30], all fail the above test. 

It should be added that the necessary condition (2.23) for having a quantum chain 
related to a solution of the Yang-Baxter equation has been examined for isotropic spin-s 
chains with s < 14 in [17, 31]. In particular, these authors found that the only spin-1 
systems satisfying (2.23) are (3 = ±1 and (3 = 00. As we have discussed above, this 
implies that for other values of /3 the hamiltonian (|4.1|) cannot be a fundamental model. 



But a priori there might exists integrable but non- fundamental models for P 7^ ±1, 00. 
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Our results, showing that there is no nontrivial three-spin charge provide a much stronger 
evidence for the non-integrabihty of all the isotropic spin-1 models with /? 7^ ±l,oo. It 
would be interesting to perform a similar analysis for s > 1; such analysis is however much 
more complicated for higher s. 

We end this section with a short remark on the general spin-s bilinear system. Is it 
possible that this system is, by a bizarre accident, an integrable fundamental model for 
some values of s? This has been answered negatively in [31] for all s < 100, using computer 
algebra. Here we present a simple calculation showing that for the bilinear systems, the 
boost operator can never produce a conserved quantity for s 7^ 1/2, which thus excludes the 
possibility of such an accident. Acting on the hamiltonian, the boost operator generates a 
candidate for the three-spin charge of the form: 



1 TT TT 1 \ ^iT^aboa ob Tiaa oh ob , ca ca r^bb oa cb r^ab i f A T Q\ 



The vanishing of this sum requires that all the terms containing spins in an arbitrary 
cluster vanish separately, which is not the case in general. The sum can vanish only if 
the terms cancel two by two, which is possible only if D"-^ ~ 5°'^, a condition which is not 
true unless s = 1/2. Therefore, the bilinear Heisenberg chain is an integrable fundamental 
model only for s = 1/2. 

5. Concluding remarks 

The simple integrability test considered in this work appears to be applicable to a 
rather general class of quantum chains. However, the range of applicability of the con- 
jectures 1 and 2 has not yet been determined rigorously. Regardless of that, this simple 
method seems to be a useful heuristic tool. It may be worthwhile to use it to try to identify 
integrable models within other physically relevant families of models. In particular, the 
completely anisotropic case of the 10-parameter family considered in section 3 remains to 
be studied in detail. However, we do not expect that it will reveal new types of integrable 
systems, beyond anisotropic generalizations of the systems found in the o(2)-symmetric 




(4.12) 



The commutator of this quantity with the hamiltonian is then 



ieA 



24 



case. In particular, we expect that the XYZ chain, the staggered XY model, the staggered 
XYZ chain (with the xy part alternating in sign), equivalent to the Dzyaloshinski-Moriya 
system (3.47), a generalized Lieb-Schultz-Mattis model, and a generalized Hubbard model 
(consisting of two copies of the XY chain interacting along their z components) are, up 
to a relabeling of variables, the only nontrivial integrable anisotropic models within the 
family (3.1). 

Conjecture 1 could be generalized in a natural way to include even models with long- 
range interactions [32], for which there seems to be no ladder operator. We would then 
require that if a hamiltonian H2 is given by a sum of two-spin interactions, there should 
exist a conserved three-spin charge H^. Observe that for models with long range inter- 
actions, the leading term in is also characteristic: although the n interacting spins in 
the leading term are no longer adjacent, the prefactor specifying the interaction of these 
n spins is distinctive. But for these models there can exist independent non-local charges 
also for finite chains (see e.g. [33]) and again it could be less obvious at first sight to assert 
that a three-spin quantity is not a product of lower order charges. The analysis of such 
models will be reported elsewhere. 

Finally, we mention a completely different integrability test for spin chains based 
on the properties of the n-magnon excitations of the ferromagnetic vacuum, which has 
been conjectured by Haldane [34]. In spin-s chains the bound-state n-magnon dispersion 
branches extend over min(Ar, 2s) Brillouin zones and integrability manifests itself in that 
all branches are real (i.e. lie outside the spin wave continuum) and are continuous through 
the zone-boundaries. This integrability criterion has been investigated for 2-magnon exci- 
tations in spin-s chains in [35, 36]. In particular, in [35] the Haldane criterion was used to 
obtain a relation between the parameters of a class of isotropic spin-s hamiltonians. This 
relation defines then a one parameter family of models, supposed to be integrable. But it 
contains only two out of the four known fundamental integrable systems for s = 3/2 [17]. 
Neither the stt(4) invariant chain [28] nor the hamiltonian Hjj of [31] belong to it. And 
apart from two points (describing the Bethe-Ansatz integrable system [37] and the chain 
related to the Temperley-Lieb algebra [38, 39, 40]), the other models within the above 
one-parameter s = 3/2 family are not finidamental integrable systems, and are very likely 
to be nonintegrable. Thus the condition in [35] turns out to be neither a necessary nor 
a sufficient condition for integrability. Presumably, further constraints should follow from 
the analysis of n-magnon excitations with n > 2. A better understanding of the Haldane 
criterion, as well as its relation to the test proposed in this work is clearly needed. Note 
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that the former approach imphes a choice of a particular vacuum, while the test is 
insensitive to the ferromagnetic or antiferromagnetic character of the model. 
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